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Physical aspects of naked singularity explosion 
- How does a naked singularity explode? - 
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The behaviors of quantum stress tensor for the scalar field on the classical background of spherical 
dust collapse is studied. In the previous works diverging flux of quantum radiation was predicted. 
We use the exact expressions in a 2D model formulated by Barve et al. Our present results show 

/<— s that the back reaction does not become important during the semiclassical phase. The appearance 

of the naked singularity would not be affected by this quantum field radiation. To predict whether 
the naked singularity explosion occurs or not we need the theory of quantum gravity. We depict 
the generation of the diverging flux inside the collapsing star. The quantum energy is gathered 
around the center positively. This would be converted to the diverging flux along the Cauchy 
horizon. The ingoing negative flux crosses the Cauchy horizon. The intensity of it is divergent only 
at the central naked singularity. This diverging negative ingoing flux is balanced with the outgoing 
positive diverging flux which propagates along the Cauchy horizon. After the replacement of the 
naked singularity to the practical high density region the instantaneous diverging radiation would 
change to more milder one with finite duration. 
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I. INTRODUCTION 

It is known that several models of gravitational collapse end in naked singularities from regular initial data. Among 
. : these models the most fully studied one is the spherically symmetric inhomogeneous dust collapse, which is described 
by the Lemaitre-Tolman-Bondi (LTB) solution [Q-g|. The final fate of this model was investigated by numerous authors 
|0-p| . They concluded that the singularity could be naked or covered by event horizon depending on the initial density 
and velocity profiles of collapsing cloud. There are also other examples of the naked singularity formation. Recent 
review of these examples and cosmic censorship hypothesis were given by some authors Pfl0| . 
S^ ' It is expected that the quantum effects analogous to Hawking radiation play an important role in the final stage of 

the naked singularity formation. From this point of view some researchers investigated a quantized massless scalar 
field on the classical background of such collapse |n|-flrj| . Their results showed that the diverging outgoing quantum 
radiation appears in the approach to the Cauchy horizon. These results are in contrast to the investigation of the 
classical radiation, i.e., gravitational waves given by the authors and K. Nakao ]l7|-p0|, where the energy flux of 
the gravitational waves remains finite as the Cauchy horizon is approached. Therefore the interpretation of this 
diverging quantum flux is important for the understanding of the nature of these naked singularities and for the proof 
or disproof of the cosmic censorship hypothesis. In another context we can also expect the intense radiation from a 
naked singularity and get information about str ongly curved spacetime. In this respect the naked singularity formation 
is discussed as a possible origin of 7-ray bursts |l5|,^l| ■ Also these observational aspects of the naked singularity might 
be useful to distinguish naked singularities from black holes on the astrophysical observations p4P3fl . 

The above quantum radiation was estimated perturbatively in the semiclassical treatment. Seemingly the divergence 
of the quantum flux would mean the breakdown of the perturbative estimate and the plausible next step will be to 
consider the backreaction of the quantum radiation. However this is not necessarily true. At the Cauchy horizon 
the classical gravity is broken down because of the causal connection to the singularity. Therefore the semiclassical 
treatment, where the gravity is treated as classical, is also broken down independently of whether or not the quantum 
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radiation diverges. It is not obvious whether the backreaction of the quantum radiation is significant or not within 
the region where the semiclassical treatment is plausible. 

There would be two criterion for deciding when the back reaction becomes important. One is that the total energy 
received at infinity becomes comparable to the mass of the collapsing star. In J24J it has been shown that, in the 
spherical dust model, the emitted flux is at most order of Planck energy before the quantum gravity phase. Here the 
quantum gravity phase means the spacetime region to where the central high curvature region can causally affect, e.g., 
the region within one Planck time up to the Cauchy horizon. If the mass of the collapsing star is much greater than 
the Planck mass, then the back reaction does not become important. The second criterion is that the energy density 
of the quantum field becomes comparable to the background energy density inside the star. To estimate the energy 
density of the quantum field we should calculate the vacuum expectation value of the stress tensor for the quantum 
field. In a four-dimensional model, it is difficult to establish this issue. However one can study the quantum stress 
tensor for a two-dimensional model which is obtained by suppressing the angular part of the 4D model |25| . For the 
self-similar dust collapse, the exact expression of it was computed |2(|. Using this model, this paper investigates the 
exact behaviors of the quantum stress tensor and discusses the generation of the diverging quantum flux inside the 
star. Also we show that the back reaction does not become significant during the semiclassical evolution inside the 
collapsing star fated to be a naked singularity. 

The plan of this paper is as follows. In Sec. [n] we introduce the calculation for the quantum stress tensor for a 



massless scalar field in the 2D self-similar spherical dust collapse performed by Barve et al. [E6|. In Sec. Ill using this 



expressions we investigate the exact behaviors of the quantum stress tensor and consider the generation of diverging 



flux and the criterion of the semiclassical treatment inside the cloud. In Sec. [V we summarize our results. We use 
the units G = c = H = 1. 

II. EXACT EXPRESSION FOR QUANTUM STRESS TENSOR IN SELF-SIMILAR DUST COLLAPSE 

We consider the exact expression for the quantum stress tensor of the massless scalar fields in a 2D collapsing 
LTB spacetime. This was first done by Barve et al. pq]. Here we briefly introduce their results, taking care of the 
consistency of signs. 

A. Self-similar spherical dust collapse 

A spherically symmetric dust collapse is represented by the LTB solution. Using the synchronous comoving coor- 
dinate system, the line element of the LTB spacetime can be expressed in the form 

ds 2 = dt 2 - A 2 (t,r)dr 2 - R 2 (t,r)dn 2 . (2.1) 

The energy-momentum tensor for the dust fluid is 

T"" = p{t,r)u»u v (2.2) 

where p(t, r) is the rest mass density and u^ is the 4- velocity of the dust fluid. 

Then the Einstein equations and the equation of motion for the dust fluid reduce to the following simple equations 

A= - R ' -. (2.3) 

, N 1 1 dF(r) 

^ B s^* (2 - 4) 

** - ^ = /W (2 - 5) 

where f[r) and F(r) are arbitrary functions of the radial coordinate r, and the overdot and prime denote partial 



derivatives with respect to t and r, respectively. From equation (2.4), F(r) is related to the Misner-Sharp mass 
function, m(r), of the dust cloud in the manner 

m(r) = 4tt / p(t, r)R 2 dR = 4tt / p(t, r)R 2 R'dr = — )-+■. (2.6) 



Hence equation ( |2.5| ) might be regarded as the energy equation per unit mass. This means that the other arbitrary 
function, f(r), is recognized as the specific energy of the dust fluid. The motion of the dust cloud is completely 
specified by the function, F(r), and the specific energy, f(r). 

The self-similar solution corresponds to the choice of 2m(r) = Xr, where A is a constant, and f(r) = 0. Following 
Barve et al. p6|, we assume that the central singularity appears at t = and rescale R as R(0, r) — (3/2) 2 ' 3 \ 1 ' 3 r at 
t = 0. The self-similar collapse solution is then 

R=(2m(r)) 1 / 3 (^(r-t)\ (2J) 



The derivative of R with respect to r is 



/9A\ 1/3 l-z/3 

fl ■ ( t) ^^ (2 - 



where z = t/r. And the first and second derivatives of R' with respect to z are given by 

dR' 2 /9A X 1/3 



dz 9 V 4 / (1 _ z ) 4 / 3 
and 

d 2 R' _ 2 /9A\ 1/3 l + z/3 



dz 2 9 V 4 / n - ^ 7 / 3 



(i- z y 



(2.9) 



(2.10) 



Using this solution, the density function p(t, r) is expressed as 



^' r ) = ^w — h* — 7S- (2 - n) 

07r(r — t)(3r — t) 



At the center the density diverges as 



Pfc0) = ^ ( 2 - 12 ) 

to the singularity. 

From the analysis of the null geodesic equations, the condition for the central singularity to be naked is obtained 

as|7} 

A< 156- 90V3 = 0.115427- ••. (2.13) 

If the singularity is naked, there exist z which satisfies 

z -^) = 2 -(t) 1/3 (i^ = - (2 - 14) 

This condition is equivalent to the existence of positive root for the following equation, 

f{v) = y i + lv i -y+\a = {) (2.15) 

where y = (1 — z) 1 ' 3 and a — (|A) . The Cauchy horizon is determined by z — z CH which satisfies this equation. 
It can be shown that this naked singularity is always globally naked jl3| . 

To obtai n th e exact expression for the quantum stress tensor, we need double-null coordinates for the 2D part of 



the metric (2.1). In the interior region, we introduce null coordinates (77, C) 

Jdz/(z-R>) (^_ i? / >0 ) 



fdz/(z-R') , D , ^ n x ' ' > 

—re J (z — R < 0) 



and 

{rJ d ^ z+R,) (z + R>>0) (2m 

^ \- r J dz ' {z+R,) (z + R'<0). { ' } 

For the convenience of the calculation for the quantum stress tensor, we introduce another null coordinates (U, V) as 

M -In £:£<$ (2 - 18) 

V-( lnC ( Z + R '>°) (219) 

V ~\ -In |C| (z + Rf <0). (Z -^> 

In these coordinates the center of the cloud is expressed by V — U from the fact that r\ = C, at the center [O. 
Hereafter we concentrate our attention to the outside of the Cauchy horizon where z — R! < 0. Then 2D metric is 
expressed as 

ds 2 = A 2 (U,V)dUdV (2.20) 

where 

,2 f[7 F1 _ f -r 2 (z 2 - i?' 2 ) (z + fl' > 0) 

In the exterior region, the metric can be expressed as 

/ 2M\ 
ds 2 = D 2 {u, v)dudv =11- — ) dudv (2.22) 

using the Eddinton-Finkelstein double null coordinates given by 

u = T-R* 1 v = T + R* (2.23) 

where 

R*=R + 2M\nf^-l) (2.24) 



and 



y/R-y/2M 



T = t- 2V2MR - 2M In —= ?-==. (2.25) 

VR + V2M 



B. Exact expression for quantum stress tensor 

For simplicity, we consider a minimally coupled scalar field as a quantum field, although the situation would not 
be changed for other massless fields, such as electromagnetic fields. 

It is known that any two dimensional spacetime is conformally flat. Then its metric can be expressed by double 
null coordinates [u, v) as 

ds 2 =C 2 (u,v)dudv. (2.26) 

If the initial quantum state is set to the vacuum state in the Minkowski spacetime ds 2 — dudv, then the expectation 
value of the stress tensor of the scalar field is given by p5| 

^>=-iM^),„ ,2 - 27) 

™~k c {h) A , < 2 - 28 > 

<T W ) = — , (2.29) 



where TZ is the two dimensional scalar curvature. 

We are now in the situation to compute the exact expression for the quantum stress tensor in the two-dimensional 
self-similar collapse. As usual, we require that the regular center is given by u = v and that v coincides with the 
standard Eddington-Finkelstein advanced time coordinate v. Assuming the relation between internal and external 
null coordinates, 



U = a{u) v = I3{V) 

we can obtain all the coordinate relations as 

v = v = /3(V) u = /3(U)=/3(a(u)). 



(2.30) 



(2.31) 



Then we can transform the quantum stress tensor given by (u, v) coordinates to the ones in the interior and the 
exterior double null coordinates. 



Using the coordinate relation equation ( 2.31 ) we can transform the equations ( 2.27 ) - ( 2.29 ) to the components 
expressed in the interior (U, V) coordinates and in the exterior (u, v) coordinates. The components of the interior 
(U, V) coordinates is given by 



(Tuu) = Fu(P') - F V {A 2 ) 
(T vv )=F v (f3')~F v (A 2 ) 

(Tuv)^-^-(lnA 2 ) 



= ±- 



247T-- " u > v 
1 z 2 -R' 2 d 2 R' 



' 24tt 2R' 



dz 2 



(z + R 1 ^ 0) 



where 



Fx{V) = ^{-^J xx 



(2.32) 
(2.33) 



(2.34) 



(2.35) 



The j3' in equation ( 2.32 ) should be considered as (3'{U). After a long calculation, we obtain 

2 



P (Ff\ l J 3 (^^ 



1 (J3'),v,v 

2 /?' 



1 



127r4(l + a;) 2 
where x is determined by the equation 

n - v 

2M 



(1 + x- 



X ,n c . X./l Z X 1 

Y } + ^ (1+ 2 } A-3^ + A + 3^ 



2 

3^3 



2 

X X 



i-21ni±^. 



(2.36) 



(2.37) 



Fu((3') can be expressed by the same equation (2.36) but with x which is determined by the equation related to the 
retarded time as 



r b - f3(a(u)) 
2M 



2 2 

+ - 



1 



3a; 3 



2 In 



1 + x 



F V {A 2 ) and F V {A 2 ) become 



1 (dR! 



1 
127T 1 4 V dz 

1 f 1 /dR' 



MA 2 ) = — {-{— l 
Fv(A 2 



12tt 14V dz 
In the exterior Schwarzschild region we obtain 



1 z 2 - R' 2 d 2 R! 

2 2R' dl? 

1 z 2 - R' 2 d 2 R' ' 



2R' 



dz 2 



(2.38) 

(2.39) 

(2.40) 



(T v . 
(T u . 



-F U (D 2 

-F V {D 2 

1 

~24tt 

24^ [? R A R 3 J 



a' 2 Fu(p') + F u (a') 



: (ln£> 2 ), 

M 2 






M\ 



The /?' in equation (2.41) should be also considered as (3'(U). 
In this region we obtain 





Fu(D 2 )=F v (D 2 ) = -—(- 1F --y 




dues are always positive for R > \M. Also we obtain 




n (n'\ 1 


a' 2 ( w 2 w\ 2 a' (w A 2w s 2w 7 w 5 \ 1 , 7 

1 III I II l°ir' 


-7u; 8 ) 


U[C1 } 127T 


4 I, 1 2A ' i) ' 4(2Af)U 'A A 3 J ' i 6 (2M ) 2 ( 


1 3w e X(-(5w 4 + 12w 3 - 8w 2 - 24w + 12)A - 6w 7 + 18w e - 12w 5 ) 





12tt 16(2M) 2 (3w 3 -2A-3w 4 -wA) 2 

where w is determined by the equation 

Th — u 2 2 1 „ , 1 — w 

- + ^r+21n . 



2M 3w 3 w w 2 
The relation between u and U is obtained as 



a'(u) 



2M 
1 — w 



2 

3w 3 



w 



w 

1 
iw 2 



(2.41) 
(2.42) 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



(2.47) 



Here w e int erpret the outgoing flux emitted from the star. For this sake, we look into the right hand side of 
equation (2.41). The first term clearly denotes the vacuum polarization, which tends to vanish so rapidly as R goes to 
infinity that there is no ne t cont ributi on to the flux at infinity. The second term is originated in the interior of the star 
as seen in equations ( 2.32J ) and ( 2.33] ) . First it appears as ingoing flux at the passage of the ingoing rays through the 
stellar surface, crosses the center and becomes outgoing flux. Since this term depends on both a 1 and /3', not only the 
outgoing but also ingoing null rays are relevant to this term. It implies that this term strongly depends on the details 
of the spacetime geometry in the interior of the star. Since the third ter m de pends only on a', only the outgoing null 
rays determine this contribution. Such a term is not seen in equation (2.32). Therefore, this contribution seems to 
originate at the passage of the outgoing rays through the stellar surface. As will be discussed later, the third term 
corresponds to the Hawking radiation, while the second term becomes important in the naked singularity explosion. 



III. EXACT BEHAVIOR OF QUANTUM STRESS TENSOR 

Using the above exact expressions, f\ we can now investigate the exact behavio rs of the q uantu m stress tensor for 
t he m assless scalar fields. The key processes are to get x which satisfy equations ( [2.37 ) and ( [2.38 ) and w in equation 

Q. ' 

From equations ( 2.23J ), (2.24) and (2.25) it is easy to confirm that the following equations hold at the dust surface, 



2M 



2M 



2 / R 

3 \2M 

2 / R 

3 V2M 



3/2 



3/2 



i R R „ > 

- l/ 2M-2M +21n 



7? 



1? 



2M 2M 



2 In 



(3.1) 
(3.2) 



Therefore x in equations( 2.37 ) and ( 2.38J ) and w in equation ( 2.4q ) become 



2M 



(3.3) 



where i? i s the value at the dust surface on the corresponding null trajectory. Therefore, to obtain the value of x in 
equation ( |2.37|) at some spacetime point, we in tegra te the ingoing null geodesic equation in the past direction from 
there to the surface of the cloud. For equation (2.38) at first we integrate the outgoing null geodesic equation in the 
past direction to the center. From there we integrate the ingoing null geodesic equation in the past direction to the 
surface of the cloud. These integration is performed numerically by the 4-th o rder R u nge-K utta meth od. We need w 
in the exterior of the dust cloud. At the surface u is obtained from equations ( 2.23 ), ( [2.24 ) and (|2.25 ). Therefore we 
should not integrate the null geodesic equation in the exterior. 

To determine what would be actually measured, the world line of the observer must be specified. For the observer 



with the velocity u^, the energy density {T^ u ) 



and energy flux (T^) 



are measured, where 



0. 



To investigate the importance of the back reaction for the central singularity formation, we compare the energy 
density observed by the comoving observer to the background energy density around the center. The energy density 
observed by the comoving observer becomes 



Pq = (T^)u»u v 



(Ti 



uu 



<r, 



VV) 



T2- 



(T, 



uvi 



r 2 (z-R') 2 r 2 (z + i?') 2 r 2 {z 2 -R' 2 ) 



(3.4) 



The results are shown in figure [j]. Basically we use the background parameters A = 0.1 and r& = fO 3 . The lines of 
p q = const are plotted in (a) and the line of p(t,r) — const is in (b). These two values coincide with each other on 
the dotted line in (b). Below this line the background density is larger than the quantum energy density. We see 
especially that near the center the background density is larger than the energy density of the scalar field until the 
central singularity. Therefore we conclude that the back reaction does not become significant during the semiclassical 
evolution. 



t The work here is based on treatment such as those in ref [E6| . One of the referees pointed out some problems with such 
an approach. In this footnote we quote his/her remarks from the ori ginal report: (i) There are reasons to believe that the 
approach is still not mathematically complete, for example, the metric (2.20) is degenerate on the CH. It is quite possible that 
the divergence observed on CH is only a manifestation of this singularity in the metric. In order for it to be taken seriously, 
one should work with an interior metric which is non-singular on the CH, relating it properly with the exterior coordinate 
system, which is not achieved so far. (ii) The coordinate invariance of the results such as divergence on CH is to be properly 
established, (iii) Approaches such as above rely on Hawking like treatment applied to the naked singular case. However, the 
situation here is significantly different when spacetime contains a naked singularity. Then the null infinity may be destroyed due 
to the presence of naked singularity. Then such a calculation has to be abandoned, especially if there are radiations escaping 
away from the stronggravity regions just prior to the epoch of naked singularity formation. In fact, calculations such as those 
of Ford and Parker |llj , Hiscock et al [[L2| are careful on this point and did consider only the marginal cases when the horizon 
and the singularity coincided, and then cutting off the spacetime. In references such as [E6|, this important point seems to have 
been ignored, and the results may be unreliable to that extent. 



W e can sh ow t hat p a does not overcome p( t, r) at the c enter for na ked case A < 0.115427 • ■ ■. Substituting equations 
(pTJ) , (U), (|j|) and (|2.10|) into equations ( |2.39|) , (|2.40|) and J2.34J ) and taking the limit r -> 0, we obtain 



F^A 2 )^^ 2 ) 



432tt 



and 



(3.5) 



<F 



c/y/ 



T 



2167T 



Since at the center v — (3(a(u)) is holds, 

Futf) = F V (P'). 
The denominators of right hand side of equation ( |3.4| ) become 

r 2 (z - R'f = r 2 (z + R'f = r 2 (z 2 - R' 2 ) = t 2 
at the center. Using these equations we obtain p q at the center as 



^ = M 2W) -2T6^ 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



The range of the i*V(/3') at the center is estimated in Appendix [A]. We conclude that p q is less than p(t, 0) at the 
center as 



Pq(«,0)< 



1 1 

< 



547rf 2 



67rf 2 



p(t,0). 



(3.10) 



One can also show that p q is not negative at the center. 

Next we consider the energy flux measured by the comoving observer. For this observer, energy flux becomes 



(T^)uV 



(Ti 



uu) 



W 



vv) 



l {z-R') 2 r 2 {z + R'f 



(3.11) 



In the interior region, outgoing flux is proportional to {Tjju) arK1 ingoing flux is proportional to (Tyv)- Here we define 
the ingoing part F; n and the outgoing part F out of the flux as 



F;, 



(T, 



vv) 



r 2 (z + R') 



A2- 



Fout = 



(Ti 



uu) 



r 2 (z-R') 2 ' 



(3.12) 



Following the sign of these two values we divide the interior region into four parts. Region 1: positive in and out flux, 
region If: negative out flux and positive in flux, region III: negative in and out flux, and region IV: positive out flux 
and negative in flux. The results are shown in figure |2|. Flux observed by a timelike observer becomes zero in region 
I and III. The observed flux vanishes at r — and dotted dashed line in figure S. We can show that the flux becomes 
zero at the center. From equations ([3Jj) and ( |3.7| ) we obtain 



(Ti 



uu) 



<T, 



vv) 



(3.13) 



Using this equation and equation (3.8) we see that the right hand side of equation ( 3.11 ) is equal to zero at the center. 
Also we see 



1 



(T uv ) = (T vv ) > — 



7 



60tt 432tt 



>0. 



(3.14) 



Therefore in and out part of the flux are positive at the regular center. 

We show this flux schematically in figure H. The intensity is proportional to the length of arrows along the r-axis. 
The comoving observers receive inward flux first and then outward flux after the crossing of the line (c). The intensity 
grows as the Cauchy horizon is approached. 

In figure [|, we show in and outgoing parts of the flux schematically. The intensity is proportional to the length of 
arrow along the r-axis. The future and past directed arrows correspond to the positive and negative flux respectively. 



The outgoing part has large intensity near the Cauchy horizon. The ingoing part has large intensity only near the 
central singularity. 

We plot these quantum values at t — const for t < in figure g and for t > in figure pi Figure [| (a) shows 
quantum energy density, energy flux, and pressure (T lil/ )n ,J "n 1 ' measured by comoving observer. The energy density 
grows positively near the center surrounded by negative energy density. The F q is negative, i.e., the inward energy 
flow and vanishes at the center. In figure |>| (b), in and outgoing parts of flux are plotted. This inward quantum 
energy flow is constructed mainly by positive ingoing and negative outgoing flux. At the center the positive outgoing 
flux cancel out the positive ingoing flux, where the net flux is zero. 

In figure || (a) we plot the quantum values at t = 0. The energy density negatively diverges as the center is 
approached. The in- and outgoing flux also negatively diverge as the center is approached. These two almost cancel 
each other. Figures @ (b), (c) and (d) denote the behaviors of quantum values for t = 0.001,0.1,1. In this region 
the energy density and the flux diverge positively as the Cauchy horizon is approached. The divergency of the flux is 
originated from the divergence of the positive outgoing flux. The ingoing flux is negative and finite. 

The same values at the surface r — r\, are plotted in figure fj]. In figure [7] (a) energy density and flux are plotted. In 
(b) in and out parts are plotted. There are inverse square divergence of the energy density, flux and out going part 
of the flux. 

We plot (Tuu), (Tyy), and {Tjjv) along the Cauchy horizon in figure g (a) and along r = r& in (b). We can confirm 
that there are positive outgoing flux along the Cauchy horizon and negative ingoing flux across it. 

Here we summarize the behavior of quantum field in the interior of the star. At first inward flux appears and then 
the positive energy is gathered near the center surrounded by slightly negative energy region. As the collapse proceeds 
the central positive energy grows and concentrates smaller region. At the naked singularity formation the gathered 
positive energy is converted to the outgoing positive diverging flux. The left negative energy goes down across the 
Cauchy horizon. It will be reach the spacelike singularity. 

Next we consider the total radiated energy received at infinity. The radiated power of quantum flux can be expressed 
as 

P = (T^u'W = cPF v {p) + F u (a') (3.15) 

where u^ is considered as the velocity of static observer. We plot this in figure g. We reproduce the results of the 
geometric optics estimates for the 4D model, i.e., the inverse square dependence for the retarded time in the approach 
to the Cauchy horizon. In this stage the behavior of P is determined by A and weakly depends on r&. On the other 
hand, during the early times, the P is proportional to the total mass M = ^Ar;, an d grows as {uq — u)~ 3 . This 
behavior agrees with the minimally coupling case for C°° background [ [l5| . 

The total energy is obtained from the integration of P with respect to the retarded time, 



(a' 2 Fu(l3') + F u (a'))du (3.16) 

at the future null infinity. The results are shown in figure |l0|. We see that the total emitted energy is much smaller 
than the mass of the original star in the range that the semiclassical approximation can be trusted. 

If we compare the a' 2 Fjj((3') and F u {a') we can say something about the origin of the diverging flux. In figure [ll] 
we plot Fu{/3') and F u (a')/a' 2 . It is seen that the first term of equation ( [3.15 ) is positive and the second term is 
negative. As was shown in p6| a' 2 behaves like 



where 



[z<mt(w) — z ut{w at the Cauchy horizon)] (3-17) 



2A 

Zout(w) = 1 - ^3 . (3.18) 



The first term of equation (2.45J) negatively diverges in the approach to the Cauchy horizon. Therefore we conclude 



that the positive divergence of the radiated power is originated from the term a' 2 Fu(f3') which has strong dependence 
on the spacetime geometry in the interior of the star. 

Let us turn to the inspection for the black hole formation. Fjj{(3') and F u (a')/a' 2 for this case are also plotted in 
figure O. For the larger u there appear the difference between the naked singularity case and the covered one. This 
difference mainly depend on the behavior of a' . In the approach to the event horizon, w goes to unity and then a' 
goes to zero. It is easy to see that only the second term of the radiated power is left in this limit and it becomes 

F u {a') = ) -J . (3.19) 

192tt(2M) 2 



This means that F u (a') is interpreted as the Hawking radiation contribution. From equation ( [2 .44 ) we obtain 



(Tw) = 7 3 (3.20) 

192tt (2M) 2 

at the event horizon. This ingoing negative flux crossing into the black hole balances with the energy loss by the 
positive Hawking flux. 

There is similar balance between in- and outgoing flux at the Cauchy horizon for the naked singularity explosion. 
In the asymptotic region positive diverging flux is observed in the approach to the Cauchy horizon. It has been shown 
above that there is negative ingoing flux crossing the Cauchy horizon. This negative flux diverge only at the central 
naked singularity. This diverging negative ingoing flux is balanced with the diverging positive outgoing flux which 
propagates along the Cauchy horizon. 

At the last of this section we give a speculation for some possible scenario of the practical quantum field radiation 
at the final stage of the collapse inspired by the above analysis. Practically the classical naked singularity should be 
replaced by an observable high-density region. After the appearance of this high-density region, the surrounding dust 
fluid falls in this region as positive energy flow. On the other hand quantum negative flux flows in this region. The 
balanced positive outgoing flux emerges from the cloud and propagates towards infinity. If the intensity of the classical 
positive energy flow is always equal to the one of the quantum negative ingoing flux, then the central high-density 
region would keep the steady state until all envelope disappears. It may be argued that in this case the inflow of 
the dust into the central region is converted to the outgoing quantum energy flux. Therefore we would observe the 
emission with finite duration rather than the instantaneous explosion. And finally the small high-density core would 
be left at the center. Even after this modification, it can be expected that the observational aspects of this radiation 
is still different from the Hawking radiation. 

IV. SUMMARY 

We have investigated the expectation value of quantum stress tensor for the massless scalar field on the 2D self- 
similar LTB spacetime. If we apply the semiclassical treatment up to the Cauchy horizon, we can depict the behavior 
of the energy flux inside the star in case of the naked singularity explosion as follows. As the dust collapse proceeds 
the quantum field flows inward and the positive energy is accumulated in the center surrounded by slightly negative 
energy. At the naked singularity formation this gathered positive energy is converted to the diverging outgoing flux. 
The negative energy envelope flows inward crossing the Cauchy horizon. The diverging outgoing flux emerges from 
the stellar surface and propagates along the Cauchy horizon. There is energy balance relation between ingoing and 
outgoing flux at the Cauchy horizon similar to the Hawking radiation. However this balance is not long-standing like 
the Hawking radiation but instantaneous. The ingoing negative flux diverges at the central naked singularity, which 
is balanced with the outgoing diverging flux along the Cauchy horizon. 

In the above picture of diverging flux, however, we omit a limitation of semiclassical treatment. The high density 
region very near the naked singularity and the part of spacetime which is causally connected to this high density region 
would not be exactly described only by the classical gravity. We need a quantum gravity theory there. Therefore 
the semiclassical treatment is not plausible there. In the suitable region for the semiclassical treatment the energy 
density of the quantum field is less than the background energy density near the center. Therefore we can conclude 
that the back reaction of the quantum radiation is insignificant to the naked singularity formation. It can be said 
that quantum gravitational effects are more important. 

In the practical point of view the naked singularity would be replaced by the sufficiently high density region (e.g., 
Planck density). In this situation the naked singularity explosion would be also replaced by a practical event. If we 
perform the replacement of the naked singularity, the inflow of the dust into the central region would be converted 
to the outgoing quantum energy flux in the assumption of the equilibrium of the classical and quantum inflows. As 
a result, we could expect that the outgoing flux would be milder and have longer time duration than the naked 
singularity explosion, in which the flux is emitted divergingly in less than a Planck time. 
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APPENDIX A: ESTIMATE OF F v (f3') 

To estimate the range of Fv((3') at the center we rewrite it as 



1 



F *W = &;\ll- W i 



2 \ 2 „2 



!(l + ») 

As we have quoted above, x is represented as 



(l+f)(ar-2) s 5 (l+f) 



3(l + x) 2 A(l+s) 



2M 



r=r b ,t=tj 



2 A 



3 (1 - Zi ) 



1/3 



(Al) 



(A2) 



where tj and Zj means the values of the junction point. For the central time t < 0, the related Zj is always negative. 
Therefore following relation holds, 



x 3 {2\ 1 2 

< 



Using the upper bound for the A we obtain 



A \3J (1-Zj) 3 



< x < 0.426. 



(A3) 



(A4) 



First we consider upper bound of the Fv(f3'). From equation (A3) it can be shown that the bracket of equation 
(Al) is less than 



fi(x)=[l- 



2 \ 2 -1 



(1 + f ) {x - 2) 2x 2 (1 + f ) 



2(1 + x) J 3(1 + x) 2 3(1 + x) 

for < x < 0.426. In this region fi(x) has a maximum value at x — as 

h(x) < /i(0) = 1. 
As a results Fy((3') is bounded as 



(A5) 



(A6) 



(A7) 



Next we consider the lower bound of Fy(/3'). The third term of the right hand side of equation ( Al) is not negative 
for < x < 0.426. We define 



f 2 (x) = 1 - 



2 \ 2 ^2 



x 2 (1 + f ) (x - 2) 



2(1 + x) J 3(1 + xY 



(A8) 



It can be easily shown that $2(x) is monotonically decreasing function for < x < 0.426. Therefore we obtain an 
inequality 



Therefore Fy((i') is bounded as 



for the naked case. 



-</ 2 (0.426) = 0.819- ■■<f 2 (x). 
5 



w. < Fv ^ * iL 



(A9) 



(A10) 
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FIG. 1. Plots of the lines p q = const in (a) and p(t,r) = const and p q = p(t,r) in (b). Each line has the following meaning. 
Solid line: singularity, long-dashed-dotted line: Cauchy horizon, dashed line: p q = p( — 1,0), dashed-dotted line: p q — 0, 
long-dashed line: p q = — p(— 1,0), long-dashed-short-dashed line: p{t,r) = p(— 1,0), dotted line: p q = p(t,r). 
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FIG. 2. The solid line and long-dashed-dotted line correspond to the singularity and the Cauchy horizon, respectively. The 
dashed line denotes the line where the ingoing part of flux vanishes. On the long dashed lines the outgoing part vanishes. 
Region I, II, III, and IV correspond to positive in and out flux, positive in and negative out flux, negative in and out flux, and 
negative in and positive out flux, respectively. On the short-dashed-dotted line the observed net flux vanishes. 
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FIG. 3. Schematic figure for the flux observed by comoving observer. Lines (a) and (b) correspond to the singularity and 
the Cauchy horizon, respectively. Line (c) denotes the line where the observed flux vanishes. 




0.5 1 1.5 2 2.5 

FIG. 4. Schematic figure for the in and outgoing parts of flux observed by comoving observer. Lines (a) and (b) correspond 
to the singularity and the Cauchy horizon, respectively, Lines (c) and (d) denote the lines where outgoing part and ingoing 
part vanish respectively. Region I, II, III, and IV are similar to the ones in figure |2| 
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FIG. 6. Plots of the quantum variables at t — 0, 0.01, 0.1, 1 in (a), (b), (c), and (d) respectively. In each subfigure upper-left 
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FIG. 7. Logarithmic plots of the absolute values of quantum variables at the surface r — rj. 
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